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ABSTRACT 

We describe a detailed analytic model for predicting statistical quantities (such as 
number counts, redshift distributions and sizes) of clusters detected in blank-field, 
thermal Sunyaev-Zel'dovich effect experiments. We include in this model the possibil- 
ity of non-Gaussian density perturbations in the early Universe and also describe a 
simple model for the effects of preheating on cluster Sunyaev-Zel'dovich effect fluxes. 
We use this model to explore the current state of the theoretical uncertainties present 
in this type of analytic modelling, highlighting where further improvement will be nec- 
essary to fully exploit forthcoming surveys. We then go on to explore the constraints 
on cosmological parameters, the presence of any non-Gaussianity and the degree of 
cluster preheating which may be obtained from both the bolocam and PLANCK ex- 
periments. We find that, providing redshifts can be measured for all detected clus- 
ters, the BOLOCAM experiment may provide detections of non-Gaussianity or preheat- 
ing and could give approximate measurements of these effects if prior knowledge of 
the various cosmological parameters is taken into account. The PLANCK experiment 
Sunyaev-Zel'dovich effect cluster survey is predicted to provide highly accurate (~ 5%) 
measurements of the degree of non-Gaussianity and preheating while also providing 
measurements of several cosmological parameters to accuracies of a few percent inde- 
pendent from those constraints that will be derived from its detections of primordial 
cosmic microwave background anisotropics. 
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1 INTRODUCTION 



The Su 



uyaev-Zelrlnvirh (R7. ) nffer-t fpurivaev K7 7el'rlnvirh 



of clusters it has the potential to reveal the effects, if any, 
of physical pro cesses such as feedback and preheating on 
the cluster gas ( iMajumdar fc Nath 2000|; ISprinKol, White & 



1972; f dr recent reviews sec Rcphacli 1995 and Birkinshaw 



Hernquist 2001; Majunidar, Nath fc Chuiba 2001) 



199£), the distortion of the cosmic microwave background 



(CMB) spectrum due to the scattering of photons from elec- 



trons ii 



hot, ionized gas (typically in clusters), is about to 



For these reasons several studies of the SZ effect have 
been made recently, using both numerical methods (e.g. | da 
Silva et al. 2000 ), analytic techniques (e.g. [Holder, Haiman 



age as a cosmological probe. To date, the SZ effect & Mohr 2001; pelabrouille, Melin fc Bartlett 2001 ) and com- 



come o: 



has b een mapped in a few clusters through targeted observa- 
tions (Carlstrom et al. 2000; Joy et al. 2001), but planned ex- 



binations of both (e.g. Kay, Liddle & Thomas 2001). Numer- 



periments such as AmiBA, ACBAR, bolocam, CBI, dasi and 
PLANCK will, for the first time, carry out statistically useful 
surveys of clusters in the SZ effect. Since the formation rate 
of clusters is strongly dependent on the values of cosmolog- 
ical parameters we may hope to place constraints on these 
parameters through such SZ surveys (Thomas fc Carlberg 



Eke, Cole fc Frcnk 1996|; |da Silva ct al. 2000|; |Rcfrcgicr ct al 



198S|; ^caramella, Cen fc Ostriker 1993t [Barbosa et al. 199'6f 



Seljak, Burwell fc Pen 2000 ; jpringcl. White fc Hernquist 



200q; [Maiumdar fc Subrahnianyan 200q; [Holder ct al. 2000 



da^ilTO~eral. 200l|; Fan fc Chuich^OOll: |Xuc fc Wu 2001 



2001| ; [Gnedin fc Jaffe 200l| ; |Kay, Liddle fc Thomas 200 iP 
Furthermore, as the SZ effect probes the baryonic content 



ical methods have the advantage of giving a full treatment 
of the complex gas physics that determines the strength of 
the SZ signal from each cluster, but suffer from being com- 
putationally expensive. Analytic methods on the other hand 
are generally highly computationally efficient, but can only 
treat the gas dynamics in an approximate way. The ideal so- 
lution would be to build a fast and accurate analytic model 
tuned to reproduce the results of numerical simulations. As 
the first blank-field SZ surveys are expected to yield results 
in the next year, the time is right to explore in more detail 
the current theoretical uncertainties in analytic SZ models in 
order to highlight where further improvement is necessary. 
It will also be crucial to ascertain just how well proposed ex- 
periments can constrain both cosmological parameters and 
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the distribution of gas within clusters. In this work we will 
address both of these questions. 

The remainder of this paper is arranged as follows. In 
§^ we describe our calculations of the thermal SZ fluxes of 
clusters. In §^ we contrast predictions for two proposed ex- 
periments, explore some of the key theoretical uncertainties 
in computing SZ survey statistics and derive the constraints 
obtainable on cosmological and gas distribution parameters 
in the near future. Finally, in §0 we give our conclusions. 



2 MODEL 



We be; in by describing how we compute the observable 



properties of SZ clusters. Our model determines only the 
thermal SZ signal, since we do not predict the peculiar veloc- 
ities of clusters needed to calculate the kinetic signal. This is 
not a significant problem since the magnitude of the kinetic 
effect is over one order of magnitude smaller than that of the 
thermal effect except at frequencies close to the null-point 



in th e the rma l effect spectral profile ( Kay, Liddle ( 
200l|). In 1 



2.1 



; Thomas 



we describe the calculation of the dark-matter 



2.2 we discuss how these halos are 



halo properties and in 
filled with baryonic material. Finally, in §2.3 we detail how 



the observable SZ properties are determined from this gas 
distribution. 



2.1 Halo Distribution 

We assume that dark-matter halos in the Universe can be 
characterised by three quantities, their mass, M, the red- 
shift at which they formed, Zf, and the redshift at which 
they are observed, Zo- From these quantities we can deter- 
mine, for given cosmological parameters, the virial radius 
and tprnperatiire of the halo from the spherical top-hat cnl- 



lapse m iodel (these properties are assumed to be fixed at the 
formation redshift and to remain fixed until the cluster is 
observed). To calculate the abundance of clusters as a func- 
tion of these three parameters we use the Press-Schechter 



theor y (Press & Schechter 1974 



Bond et al. 



1991 ), plus exten sions described by Lacey & Cole (199J) 



199l|; ^ 



and Sasaki (199^ ) that allow the distribution of formation 
redshifts to be determined. 

The Press-Schechter theory predicts that the number of 
clusters of mass M to M -|- dM per unit comoving volume, 
V, observed at redshift Zo is 



dMdV 



/Pcrit^o dy 
M dM 



P{y), 



(1) 



where pcrit is the critical density of the Universe at the 
present day, y — 5c{z)/a{M), with 5c{z) being the critical 
linear overdensity for collapse at redshift z divided by the 
growth factor (which is normalised to unity at the present 
day), a{M) being the linear theory mass variance at the 
present day in a sphere containing mass M on average, 
P{y) being the probability distribution function for density 
perturbations in the early universe and f — 1/ P{y)dy. 
Typically, P{y) is taken to be a Gaussian (as expected from 
traditional models of infiation), but we will also consider 
other forms in this paper. This distribution is easily ex- 
pressed in terms of Zc, using the volume-redshift relation for 
Friedmann cosmologies 



^=Anr'{zo)cH-\z,), (2) 
dzo 

where c is the speed of light, r{z) is the coordinate dis- 
tance to redshift z, H{z) is the Hubble constant at that 
redshift and we have computed the volume over the whole 
sky. Therefore 



d^n 



d^n dV 



dMdzo dMdV dzo ' 

Finally, we consider two models for the distributi on of halo 
formation redshifts. In the model of ^asaki (199^ ) the frac- 
tion of clusters of mass M seen at redshift Zo which formed 



2001) 



betw een redshifts Zf and z{ + dzc is given by (Verde et al 
[ dP(yf) dy 



dz{ 



P(yo) dy dzf ' 



(4) 



Therefore, the number of clusters seen in the range M to 
M + dM, Zf to Zf + dZo and z^ to Zf + dzf is 



d-^n 



d^n dV df 



(5) 



dMdzodzc dMdV dZo dzi 

An alternative derivation of the distribution of halo for- 



mation times due to Lacey & Cole (1993) gives the result 

yt 



dzi 



Jm/2 

(j{M) 



1/2 M 



M/2 

2 dlna 



M' {a{M'Y/G{MY - 1)3/2 



cr(M')2 dluM 



X exp 



(yf - i/o)V2 

cr(M')2/cr(M)2 - 



dM', 



(6) 



for the case of Gaussian initial conditions. The difference 
between the two results is due to the different way in which 
"forrn ation" is defined in the two approaches. Lacey & Cole 



(1993| ) define the formation event as that time at which a 



halo first had a progenitor at least half as massive as itself. 



Sasaki (1994) instead derives expressions for the formation 



and destruction rates of halos under the assumption that the 
destruction rate has no preferred mass scale. Figure |l] com- 
pares the two formation-redshift distributions for clusters of 
masses lO^*, lO"'^ and lOi^/i'^M© (thin, heavy and very- 
heavy lines respectively) observed at Zo = in a ACDM 
cosmological model. Note that although the shapes of the 
two distributions differ, both decline rapidly at high red- 
shifts and the mean redshift of formation becomes lower as 
the cluster mass (and hence j/o) increases. 

It is apparen t that the two distrib utions are rather dif- 
ferent, with the l^^acey fc Cole (1993 ) model predicting a 
much higher mean redshift of formation. The question of 
which is th e best distribution is almost certainly problem 
dependent. Lacey & Cole (1994) show that their calculation 
accurately describes the results of N-body simulations given 
their defini tion of forma tion time. Note that the expression 
derived by Sasaki (1994[) is q ualitatively similar to that de- 
rived by Lacey & Cole (1993) using t he mass-halving time s 
for single trajectories, although the Lacey & Cole (1993) 



expression contains a dependence on the shape of the dark- 
matter power spectrum whereas that of ISasaki (199^ ) does 
not. However, for our work we are really interested in know- 
ing the properties of gas in a halo. Halos of given mass that 
formed earlier will be denser and hotter. In reality, the pro- 
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Figure 1. Comparison of formation redshift distributions. The 
curves show the fractional formation rate as a function of red- 
shift for a clusters of masses 10^'', lO^'' '^ and lO^^h'^MQ (thin, 
heavy and very heavy lines respectively) observed at Zo = in 
a ACDM cosmology. The corresponding value of j/o is given for 



each cluster in the figure. The solid lines indicate the Laccv 
Cole pi 9^ calculation while the dashed lines indicate the ^asa ki 
(1994|)calculation. 



cess of cluster formation is a continuous one and so it is 
not clear just how the formation redshift should affect the 
gas properties. In this work we will therefore consider both 
distributions, and also the case of Zf = Zo. Two of these 
possibilities are limiting cases. The case Z{ — Zo is clearly 
the lower lim it on formation redsh ift, while the distribution 
proposed by Lacey & Cole (1993) should give a reasonable 
upper limit as it is difficult to see how the observed gas 
properties could be determined by the physical state at any 
earlier time (i.e. when less than half of the final mass of 
the cluster is in place). Our default distribution will be the 
Sasaki model. 

With P{y) a Gaussian the Press-Schechter theory pre- 
dictions differ significantly from the results of numerical sim- 
ulations of structure formation (by a factor close to 2 at the 
characteristic mass Af*). Several fitting formulae have been 
proposed which produce a much better agreement with nu- 
merical results. We will consider two such fittin g formulae, 
those proposed by Sheth & Tormen (1999) and Jenkins et 
al. (200^ (hereafter ST and J2000 respectively), which can 
be described by 



Pst{v) = 0.3223^/^ 



1 + 



(ay 



2\q 



exp 



-Pj2ooo (y) = 



A 



■ exp I 



In 



y 



1.686 



+ c 



) (7) 



(8) 



where a = 0.707 and q — 0.3 for the ST formula and A = 
0.301(307), c = 0.64(0.61) and p = 3.88(3.82) for the J2000 
formula in the ACDM(tCDM) cosmology considered in 

To explore non- Gaussian initial conditions we make 



use of the log-normal probabilit y distribution proposed by 
Robinson, Gawiser fc Silk (2000|), 



PRGsiy) = 



fC 



2tt\A 



exp 



- l^Hy) 



(9) 



where x(y) = \n(B + Cy\A\/A) /\A\, with yl, B and C as 



defined by Robinson, Gawiser & Silk (2000). The degree of 



non-Gaussianity is fully specified by the parameter A (with 
A — corresponding to the Gaussian limit). We will, how- 
ever, characterize non-Ga ussian models by the more intu - 
itive variable G, defined by Robinson, Gawiser fc Silk (200C ) 
as the number of > 3cr peaks in the non-Gaussian distribu- 
tion relative to the number in a Gaussian distribution (such 
that G — 1 represents the Gaussian limit). 

2.1.1 Cosmic Variance 

The above calculations allow us to calculate the mean abun- 
dance of dark-matter halos as a function of their mass, and 
redshifts of observation and formation for any given cosmo- 
logical parameters. However, for SZ surveys covering rela- 
tively small fields of view it is important to assess the effects 
of sample variance (a.k.a. cosmic variance) in order to accu- 
rately determine the ability of the survey to discriminate be- 
tween models. In fact, the sample variance is expected to be 
rather small (close to Poissonian) since the SZ effect probes 
a wide range of redshifts such that any intrinsic correlations 
between clusters are diluted (as can be estimated using an- 
alytical calculations of the cluster bias to compute their an- 
gular correlation function). To measure the sample variance 
we make use of the Hubble Volume simulations which were 
carried out by th e VIRGO Consortiu m and which are pub- 
lically available ( Evrard et al. 1998 ). These large N-body 
simulations have been used to construct catalogues of dark- 
matter halos (listing the mass and observed redshift) as seen 
along a past lightcone. 

To compute the cosmic variance from these simulations 
we construct mock surveys of the required angular size by 
choosing a random line of sight through the simulation light- 
cone and assigning each cluster in the field of view the SZ 
flux computed from our model for a cluster of the same mass 
and observed redshift and with a formation redshift drawn 
at random from the distributions described above (note that 
this assumes that the spatial distribution of halos of given 
mass and observed redshift is independent of their forma- 
tion redshift, which is unlikely to be correct in detail). We 
then compute the statistic of interest from this mock survey. 
Repeating for several mock surveys we obtain a measure of 
the variance in the statistic. Results from these calculations 
will be shown in 



3.3.1 



2.2 Baryon Distribution 

2.2.1 No Preheating/No Cooling 

Numerical simulations of gas in CDM clusters, for exam- 
ple th^SaiitaBajbara cluster comparison project (SBCCP) 
of Frenk et al. (1999), show that, in the absence of heat- 



ing and cooling processes (other than shocks and adiabatic 
compression) the baryonic material of a cluster has a density 
profile well de scribed by a beta-model (Cavaliere & Fusco- 
Femiano 1976 ), and is close to hydrostatic equilibrium, with 
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a temperature at the virial radius equal to approximately 
half the virial temperature of the cluster. For each cluster in 
our calculations we determine the virial temperature from 
the spherical top-hat collapse model and the total mass of 
gas using Mgas = {ilh/^o)M . This gas is distributed in the 
cluster as a beta-model with /3 — 0.9 and a core radius equal 
to 0.09 of the virial radius. We set the gas temperature at 
the virial radius to one half of the virial temperature and 
determine the temperature at smaller radii by assuming hy- 



L equihbimm. As 



TTS- ell the 



drostat: 

profile bf the SBCCP cluatcr. Beyond the virial radiua wc 
extrapolate the pressure and density as power-laws which 
also provides a good fit to the SBCCP cluster (hydrostatic 
equilibrium is not a good assumption here as infall is occur- 
ring). 



2.2.2 Preheated and/or Cooled Profiles 

It is well known that, in the absence of radiative cooling 
or any other process (other than infall-driven shocks) the 
relation between the X-ray luminosity and temperature of 
clusters differs significantly from that observed. One possi- 
ble solution is that the intra-cluster medium (ICM) gas was 
preheated before it collapsed into the cluster, increasing its 
entropy and resulting in a shallower density pro file. (For fur - 



ther discussion of these points see, for example, Fozzi 2001 



In the case of no preheating the entropy of the cluster 
gas arises almost entirely from the shock it experiences as 
it accretes onto the cluster. If we assume that the gas is 
in hydrostatic equilibrium, then its "entropy" density as a 
function of radial distance r from the cluster centre is given 
by 



Kshock{r) 



1 



GM{r')p,{r') ^^, 



+ 



Pg(' vir, 

^niH 



■kBT'g(rvir) 



(10) 



where pg and Tg are the gas density and temperature re- 
spectively, 7 = 5/3 is the ratio of specific heats for the gas, 
M{r) is the cluster mass profile, rvir the cluster virial ra- 
dius and fi = 0.59 the mean atomic mass appropriate to a 
fully ionized primordial gas (with hydrogen mass fraction 
X = 0.76). 

Motivated by X-ray studies (Ponman, Cannon & 



Navarro 1999) we introduce a minimum entropy such that 
the entropy profile becomes. 



K{r) 



-^min ^ *^ ^crit 

-R'shock(r) r > r-crit 



(11) 



where rciit is defined by A'shock(rcrit) = -K'min. Assuming 
hydrostatic equilibrium and an NFW density profile for both 
dark matter and shock-heated gas the corresponding density 
profile of the preheated gas can be obtained (by expressing 
the equation of hydrostatic equilibrium in terms of pg and 
K, and then solving for pg) and is given by (S.-P. Oh, private 
communication) 



(§^PNFw(r)) 

X ^ln(l + CNFW) 

^ n^PNFw(r) 



7 — 1 Vs 

^NFW 

(1 + CNFw) 

2 _ ln(l + x) ^ 



(12) 



r < rcrit 
r > rcrit 



where cnfw is the concentration parameter of the NFW 
profile (defined as the ratio of viri al to scale radii and com 
pute d using the method given by 



Navarro, Frenk & White 



1997). From this the corresponding pressure and tempera- 



ture profiles are easily found. The effect of the minimum 
entropy is to reduce the central density and pressure of 
gas in halos, the effect being strongest in lower-mass ha- 
los. This necessarily reduces the SZ fiux of the halos and so 
affects the abundances and redshift distributions of SZ se- 
lected samples. We choose to express the minimum entropy 
density as K^i n = 10^**A'34/i"*^^ergs cm"^ g"^''^. [Tozzi & 



Norman (200l| ) found that values of K34 ~ 0.12 and 0.25 



were required to explain the properties of groups and clus- 
ters with X-ray temperatures below and above 2 keV respec- 
tively, or alternatively that a redshift dependent value of 

— 0.5/(1 -I- z) fitted the data over the whole tempera- 
ture range for which data was available^ 

There remains some doubt as to whether preheating 
is necessary at all. Dave et al. (2001) use hydrodynamical 
simulations which include radiative cooling to show that a 
good fit to the cluster X-ray luminosity-temperature relation 
can be achieved when cooli ng is included without the need 
for any preheating (see also Pearce et al. 2000; Bryan 200C; 
Muanwong et al. 200ll) . 



Unfortunately, the effects of cooling on cluster gas pro- 
files are more difficult to ascertain analytically, and this will 
be a key area in which analytic modelling must improve 
to fully exploit future SZ surveys to constrain the phys- 
ical processes affecting cluster gas. The high-density /low- 
entropy inner regions of the cluster should cool in times 
much shorter than the Ifubble time. Current evidence from 



numerical simulations (e.g. Pearce et al. 2000; Muanwong et 



al. 2001) suggests that as gas cools the surrounding higher 



entropy gas ffows inwards to take its place. If this is indeed 
the case the cooling will also introduce a minimum entropy 
into cluster gas profiles. As such, we expect our model of pre- 
heating to produce qualitatively similar hot gas profiles to 
those resulting from radiative cooling. Therefore, lacking at 
present any better way to parameterise the effects of cool- 
ing we will consider preheated models only, and recognise 
that what looks like preheating may in fact be the result of 
radiative cooling. 



2.3 SZ Calculation 



Our calculation of the SZ flux essentially follows that of Kay, 



Liddle & Thomas (2001). Having determined the density 



and temperature proflles of the electrons in each cluster we 
determine the Compton parameter as a function of angular 
radius, y{8), using 



^ mm 



10'^'*K34 ergs 



1 Note that [Tozzi fc Norman (200l[ ) define 

cm~^ g~^^^, i.e. without the explicit h dependence of our defini- 
tion. We have adjusted their values of K34 accordingly. 
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Table 1. Parameters of the two experiments, BOLOCAM and PLANCK, considered in this work. 

Experiment Solid angle (sq. deg.) Beam FWHM Central frequency Frequency response Background, j/j^g + Scrj^g 

BOLOCAM 1 1' 143 GHz Top-hat, Au/u = 0.07 1.8 x 10"^ 

PLANCK 41253 (all sky) 8' 143 GHz Gaussian, a = 22.5GHz 1.0 X IQ-^ 



moC^ 



(13) 



where fee is Boltzmann's constant, ctt is the Thomson 
cross section, mo is the mass of an electron and the inte- 
gral is taken along the line of sight passing a distance 9 
away from the centre of the cluster. The electron density 
is ric — 0.88pg/mH, where mn is the mass of a hydrogen 
atom, as appropriate for a fully ionized primordial gas (with 
hydrogen mass fraction X — 0.76). Although the above inte- 
gral extends from — oo to -|-oo, in practice it converges with 
limits of order the virial radius of the halo. This profile is 
then convolved with a beam to mimic instrumental effects, 
giving 



y{e')bi\e-e'\)d'e', 



(14) 



where b{d) is the beam profile, which we take to be a Gaus- 
sian. The total SZ flux increment or decrement at frequency 
v is then found using 



= 2-rSoG{xo) 



ysio)ede, 



(15) 



where 5*0 = 2.29 x 10'* mjy arcmin"^ with 9 in arcmin and 
9h is the radius at which the Compton y parameter of the 
observed cluster profile falls below the background. We de- 
fine 

/>CXD 

G{xo) = / g{x)f{x - xo)dx, (16) 
Jo 

where x = hiz/kBTcMB, h is Planck's const ant, Tcmb 



2.72 5K is the mean temperature of the CMB (Mather et al 
19990 



e" + 1 



(17) 



is the frequency dependence of the thermal SZ effect and 
f{x) is the filter response function, which we will take to be 
either a Gaussian or a top-hat. 

2. 3. 1 Background 

The limiting flux for SZ detections will depend upon the SZ 
background due to the summed contributions of all unre- 
solved sources. This background can be characterized by a 
mean Compton y parameter and an rms fluctuation about 
this value. We can calculate the background and its fluc- 
tuations fro m our model using a n approach similar to that 
described by Bartelmann (2001) (similar calculations have 



been used to compute the power spectrum of the SZ effect 



by Bruscoh et al. 2000; Valageas, Balbi fc Silk 2001 



Benson 



et al. 2001 ). We compute the smoothed distribution, ys{9), 
for each point in the (M, Zo,Zf) plane and then determine 
the mean j/s, y, and its variance, a^, for each cluster, av- 
eraged over the whole sky. Assuming the positions of the 



clusters on the sky to be uncorrelated (a good approxima- 
tion as discussed in ^ 2.1.1), the central limit theorem implies 
that, when the contributions from all clusters are added to- 
gether, the observed mean background will be i/bg = 
with variance (Tbg = of, where the sums are taken over 
all clusters^. 

We take the background to be i/bg + 3(Tbg since below 
this any signal will be lost in noise due to the unresolved 
halos. The computed values of the background will differ 
between experiments, since the value of abg depends upon 
the beam width, and also depend upon the cosmological 
model. (Note that our values of j/bg are consistent w ith the 
upper limit of 1.5 x 10~^ from the firas instrument; 



et al. 1996 



Fixsen 



,) Typically we find backgrounds in the rCDM 
cosmology to be approximately half th ose in the low-densit y 
models (comparable to the results of |da Silva et al. 2000 ). 
Since the effective background will depend to some degree on 
the particular analysis method applied to data we choose to 
adopt a single value of the background for each experiment 
(i.e. the background is assumed independent of cosmology) 
to permit easier comparison between models. These values 
are listed in Table |l|. 



3 RESULTS 

We pre sent results tuned specifically to the BOLOCAM exper- 
iment ( Glenn et al. 199"8[) but will also consider the PLANCK 
experiment (Passvogel & Felici 200C) which will survey a 



much larger area, albeit with lower angular resolution than 
BOLOCAM. For our purposes, these two experiments differ 
mostly in terms of the solid angle of sky which they will 
survey and their beam size which determines the flux of 
each cluster, and also influences the background. The spe- 
ciflc parameters used for these two experiments are listed in 
Table |l|. 

We will also consider three different cosmological mod- 
els, ACDM (which is currently the most favoured obser- 
vationally), rCDM and OCDM. The default parameters 
of these models are given in Table ^ and are used for 
all calculations unless otherwise noted. We choose to flx 
the value of Q,hh? (we deflne the Hubble constant to be 
Hp = lOOfakm/s/Mpc ) to 0.019 based on the results of 
Buries & Tytler (199J). The val ue of as is flxcd from th e 



observed abundance of clusters (Eke, Cole & Frenk 1996), 
and the spectral shape parameter, F, is chosen to match the 



shape o f galaxy power spectra on large scales (Percival et 
al. 2001). All the models have a primordial power spectrum 
with spectral index n = 1. 

We now present results for key observable properties of 

^ Since the distribution of y is different for each cluster the cen- 
tral limit theorem only ap plies under certain conditions (see, for 



example, Eadie et al. 1971), but we have checked that these apply 
for this particular calculation. 
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BOLOCAM = 
PLANCK : 



|S„5t„J>l mJy 

BOLOCAM 
PLANCK 



0-1 - 
I 
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ISJ/mJy 



§ 1 ^ 



|S„so„J>l mJy , 
BOLOCAM : 
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Figure 2. Properties of SZ clusters as detected by BOLOCAM (dotted lines) and PLANCK (dashed lines). The left-hand panel shows counts 
as a function of flux at 143GHz, the centre panel shows the redshift distribution of clusters brighter than | Si 43GHz I = ImJyi and the 
right hand panel shows the distribution of angular sizes for the same clusters. All results are for the ACDM cosmology with the J2000 

P{y)- 



Table 2. Parameters of the cosmological models considered in 
this work. 



Model 


no 


Ao 




h 


0-8 


r 


ACDM 


0.3 


0.7 


0.0388 


0.7 


0.93 


0.21 


rCDM 


1.0 


0.0 


0.0760 


0.5 


0.52 


0.21 


OCDM 


0.3 


0.0 


0.0388 


0.7 


0.87 


0.21 



SZ clusters. We begin by examining the dif fere nces between 
the BOLOCAM and PLANCK experiments (§3.1), proce ed to 
explor e the t heoretical uncertainties in the results (§p.2.l|, 



3.2.2 



3.2.3), exam ine b riefly the dependence on cosmo- 
3.3.1) and finally consid er th e effects 



logical parameters 



of preheating (§3.3.3) and non-Gaussianity (§3.3.2), focus- 
ing in particular on what constraints can be put on these 
mechanisms and also on cosmological parameters. 



3.1 A Case Study: BOLOCAM vs. PLANCK 

We begin by briefly contrasting the two experiments that 
will be considered in this work. (We show results for the 
ACDM cosmology and the J2000 P{y).) Figure | shows 
three statistics of the SZ clusters that may be seen by 
the BOLOCAM and PLANCK experiments. The left panel is 
a plot of the number of clusters expected versus the flux 
limit. The middle and right panels are plots of the redshift 
and angular size distributions respectively for clusters with 
|S'i43GHz| > ImJy. The redshift distribution is normalized 
by the total number of clusters brighter than this flux. (Un- 
less noted otherwise, the format of all subsequent figures 
showing SZ survey statistics will be the same.) 

It is immediately obvious that the smaller beam of the 
BOLOCAM experiment allows it to detect many low-flux clus- 
ters which would disappear into the background for PLANCK. 
As a result, bolocam should flnd over ten times more cluster 
per square degree than PLANCK, resulting in approximately 
10 SZ cluster detections per square degree. Of course, the 
full-sky coverage of PLANCK will mean it detects many more 
clusters (several thousands) in total. Because of its greater 
angular resolution, BOLOCAM will detect clusters with a 



much broader redshift distribution, which will help discrimi 
nate betw een models with different cosmological parameters 
(see 



3.3) 



Finally, the distribution of cluster sizes also re- 
flects the size of the experimental beam in each case, with 
only a small fraction (approximately 10% for BOLOCAM and 
1% for planck) of detected clusters being resolved (as was 
noted for the planck experiment by Kay, Liddle & Thomas 



2001) 



Given the predicted slope of the SZ cluster counts we 
find that the currently proposed BOLOCAM survey (1 square 
degree surveyed to a limit of ImJy) is close to optimal for 
fixed total observation time in terms of maximizing the num- 
ber of detections. A slightly shallower and broader survey 
would detect somewhat more SZ clusters, but a deeper and 
narrower survey has little advantage since the cluster counts 
begin to turn over because the small, faint clusters begin to 
blend into the background due to the beam size. 

Figure ^ shows which regions of the Afhaio-Zf plane are 
detectable at different observed redshifts for a | Si 43GHz | ~ 
ImJy limited survey. Each line in the figure corresponds to 
a single Zobs, and only clusters with (Mhaio,Zf) above and 
to the right of the line would be detected. It is immediately 
apparent that the minimum cluster mass detectable depends 
only weakly on the observed redshift. The dependence arises 
through the angular diameter distance (since the total SZ 
fiux scales approximately as d^^ in the absence of instru- 
mental beam effects). The minimum mass for detection is 
also seen to be lower for objects with higher since these 
objects are denser and so produce larger SZ flux for given 
mass. (The lines become horizontal at Zf — Zohs since no 
cluster can have formed after it was observed.) 



3.2 Assessing the Theoretical Uncertainties 

In this subsection we explore the theoretical uncertainties 
inherent in current analytic models of SZ cluster properties. 
Our aim is to assess the relative importance of these un- 
certainties and highlight where further work is needed to 
produce an accurate model. 
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Figure 3. The regions of tlie Mjjaio^Zf plane in wliicli clusters 
can be detected in an |S'i43ghzI = ImJy limited BOLOCAM survey. 
For clusters observed at given 2:obs (as indicated by the different 
line types), only those with (Mhaio) ^() above and to the right of 
the line are detected. 



3.2.1 Sensitivity to the Mass Function 

Many previous analytic studies of the SZ effect have use d the 
Press-Schechter mass function (Press & Schechter 1974), but 



more accurate fitting formulaes are now known. To explore 
the differences intriduced by these alternative formulae, we 
have implemented both the Press-Schechter mass function 
and two others based upon N-body simulations developed 
by J2000 and ST. The J2000 mass function is currently the 
best match to numerical simulations over a wide range of 
mass scales, and should be taken as giving the best estimate 
of the halo mass function. We show results for the other 
mass functions simply for comparison. 

The results for all three mass functions for the BOLO- 
CAM experiment are displayed in Fig. ^ The main differences 
between the models show up in the redshift distribution of 
clusters. The ST mass function produces the largest frac- 
tion of high redshift clusters while Press-Schechter produces 
the least (a factor of 4 fewer than ST for z > 1). Changing 
the mass function also produces slight vari ation s in the total 
number of clusters seen. As we will show (§3.3), an accurate 



redshift distribution will be crucial in determining cosmolog- 
ical parameters and details of the gas physics, in particular 
for the BOLOCAM experiment. For that experiment, number 
counts alone do not contain sufficient information to be able 
to provide interesting constraints. 

The differences between models with these three forms 
for P{y) as shown in Fig. ^ should be compared to the 
expected statistical errors from individual experiments. In 
the case of BOLOCAM the statist ical errors (which should be 
close to Poissonian — see §3.3.1) will be significantly larger 
than the uncertainties due to the form of P{y). For the 
PLANCK experiment however, the statistical errors will be 
much smaller, making the systematic differences due to the 
form of P(y) the dominant uncertainty. 



3.2.2 Dependence on Formation Redshift Prescription 

We next examine the effect of varying the prescription for 
cluster formation redshifts. The formation redshift deter- 
mines the cluster's gas density and pressure, which directly 
impact the detectability of the cluster. This is a potentially 
significant effect on the number counts of SZ clusters. 

At present it is not fully clear how a cluster's formation 
and merger histories influence the thermodynamic proper- 
ties of the gas it contains. While numerical simulations may 
soon clarify this issue we have for now explored three pos- 
sibilities for the formation redshifts of clusters (which we 
assume fix the thermodynamic properties of the gas as de- 
scribed in ^). 

Comparing the predictions of the three formation red- 
shift forr nulas in Fig. we find that the fun ctions pro- 
posed by Lacey & Cole (1993) and Sasaki (1994) yield quite 
similar results. However, the lower limit of Zf — Zo yields 
substantially smaller number counts. Interestingly, all three 
models produce remarkably similar normalised redshift dis- 
tributions. It is worth noting that the differences in num- 
ber counts for our three formation-redshift distributions are 
greater than the differe nces r esulting from the various forms 
for P{y) discussed in j j3.2.1 , and at ImJy become compa- 
rable to the random errors expected in the BOLOCAM ex- 
periment. These differences highlight the need for a better 
understanding of how a cluster's gas properties are deter- 
mined in order to make sufficiently accurate calculations of 
the abundance of SZ clusters, and to allow cosmological and 
gas-distribution parameters to be determined without sys- 
tematic biases. 



3.2.3 Dependence on Gas Temperature Profile 

The form of the assumed gas density and temperature pro- 
files will affect the calculation of the SZ flux — here we ex- 



amine three te mperature profiles. (Note that Kay, Liddle & 
Thomas (2001) explored the effect of changing the gas den- 
sity profile, finding that a beta-model profile as used in this 
work produces SZ fluxes around 40% higher than an NFW 
gas profile.) In our standard case, the outer boundary tem- 
perature was set by the virial theorem, the density set to a 
beta-model, and the cluster's temperature profile calculated 
assuming hydrostatic equilibrium (as described in §2.2.1). 
Since the resu lting mass-temperature relation disag rees with 
observations (Horner, Mushotzky & Scharf 1999), we also 



considered a model in which the observed mass-temperature 
relation was used to determine the outer boundary temper- 
ature with the temperature profile again calculated from hy- 
dr ostatic equilibrium. Specificall y, we used the relation given 
by Kay, Liddle fc Thomas ( |2001 ) which was based upon th e 



observational data of Horner, Mushotzky & Scharf (1999) 



Finally, the sensitivity of the SZ flux to the temperature 
profile was tested by using an isothermal profile. While this 
does not result in hydrostatic equilibrium, it should provide 
a useful check as to the sensitivity of the SZ flux to the tem- 
perature profile. In this case the temperature of the gas was 
set equal to the virial temperature of the cluster. 

As seen in Fig. ^ the effects of these different temper- 
ature profiles are very small, about 20% (i.e. significantly 
smaller than the uncertainties due to P{y) or the forma- 
tion redshift distribution, and less than the 40% differences 
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Figure 4. Properties of SZ clusters for different determinations of the halo mass function, as characterised by P{y). Solid, dotted and 
dashed lines show results for the PS, ST and J2000 forms for P(y) respectively. The left-hand panel shows counts as a function of flux at 
143GHz, the centre panel shows the redshift distribution of clusters brighter than |S'i43ghzI = Inijy, and the right hand panel shows the 
distribution of angular sizes for the same clusters. (All subsequent figures will follow this general format.) All results are for the ACDM 
cosmology and the parameters of the BOLOCAM experiment. 
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Figure 5. Properties of SZ clu sters for different formation redshift calculations. Dotted and dashed lines show results for the Sasaki 
(1994Qnd ^ acey fc Cole (I993I ) formation redshift distributions respectively, while the dot-dashed line indicates the results when halos 
are assumed to have formed at the redshift they are observed at. The left-hand panel shows counts as a function of flux at 143GHz, the 
centre panel shows the redshift distribution of clusters brighter than |Si43GHz| = Inijy, and the right hand panel shows the distribution 
of angular sizes for the same clusters. All results are for the ACDM cosmology and the parameters of the BOLOCAM experiment and use 
the J2000 P{y). 



found by Kay, Liddle & Thomas (2001) for different gas 



density profiles), a consequence of the fact that the temper- 
ature varies only slightly with radius even in the hydrostatic 
equilibrium profiles. Nevertheless, improvement in this area 
would be necessary to fully exploit the high statistical accu- 
racy of SZ cluster counts obtainable by the PLANCK experi- 
ment. 



3.3 Constraints on Cosmological Parameters, 
Non-Gaussianity and Preheating 

3.3.1 Dependence on Cosmological Parameters 

It is well known that the abundance and redshift distribution 
of SZ clusters is sensitive to cosmological parameters as was 
discussed in our Introduction. Here we will briefly examine 
results for different cosmological parameters to demonstrate 



how SZ survey statistics are affected by these quantities , and 
then explo re the constraints that may be obtained in §3.3.2 
and 53.3.3. 



We consider three different sets of cosmological param- 
eters (defined in Table ||), rCDM, ACDM, and, for com- 
pleteness, OCDM which is now strongly rule d out by CMB 
measurements (e.g. de Bernardis et al. 2001). Fig. M shows 
that there are large differences in both the total number 
counts of clusters and in the redshift distribution between 
rCDM and ACDM. As noted before, the large differences 
in the redshift distributions imply that obtaining redshifts 
for SZ detected clusters will be very valuable in helping to 
constrain cosmological parameters. The errorbars in the left- 
hand panel of Fig. [7| show the rms scatter obtained from our 
cosmic variance calculations. As expected, the scatter is very 
close to Poissonian since the wide SZ redshift distribution 
projects out most of the clustering signal. 
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Figure 6. Properties of SZ clusters for different assumptions about the temperature profile of cluster gas. Solid, dashed and dotted lines 
indicate our standard assumptions, isothermal halos and halos given a temperature from the observed M-T relation respectively. The 
left-hand panel shows counts as a function of flux at 143GHz, the centre panel shows the redshift distribution of clusters brighter than 
|Si43GHzl = ImJyj and the right hand panel shows the distribution of angular sizes for the same clusters. All results are for the ACDM 
cosmology and the parameters of the BOLOCAM experiment and use the J2000 P{y)- 



The most significant difference is in the redshift distri- 
bution. Even a few clusters at z > 1 will rule out f2o = 1 in 
the absence of significant non-Gaussianity, providing some 
prior constraint on erg is adopted (since even with redshifts 
t he co nstraints on f^o and erg are highly degenerate, see 
§3.3.2). The results are also sensitive to erg. Varying erg 



slightly scales the number counts up and down significantly, 
and alters the high-z tail of the redshift distribution. 

The observed SZ flux will also depend on the value of 
Sib- In our model, where the effects of radiative cooling are 
ignored, the electron density (and hence y{6) and ys{0)) in 
each dark matter halo scales in proportion to fib. Further- 
more, the background which our model predicts also scales 
in proportion to fib- As a result of these two scalings ^b 
in eqn. ( |l5| ) is independent of f2b and so the observed, inte- 
grated SZ flux of a model cluster simply scales in proportion 
to rib- 



3.3.2 Constraints on Non-Gaussianity 

As shown in Fig. ^, non-Gaussian initial conditions will have 
a large effect on both dN/dz and A''(> S). Cluster sizes are 
mostly unaffected. These effects are degenerate with changes 
due to varying several of the cosmological parameters. We 
will now explore quantitatively these degeneracies, and the 
constraints that can be placed on G and various cosmological 
parameters. 



Table 3. The priors used to obtain constraints on G and K34. For 
each parameter we assume a prior probability which is a Gaussian 
with a equal to the fraction given in the table. 

Parameter "True" value Prior accuracy 



Ho 

Hb 
o-g 



0.3 
0.04 
0.9 



15% 
20% 



not necessarily the current best-fit values, for this analy- 
sis). For flo and fib we adopt priors based on the analysis 
of CMB anisotropics of (Netterfield et al. 2001, specifically 
their "Flat, LSS & SNla" determinations), while for erg we 
use the fractional errors obtain from cluster abundance anal- 



ysis by Eke, Cole , 



Frenk (1996). We will show results both 

model with 



with and without these priors. We will take 
G 



(2OO0I ) 



5 as an example (note that Robinson, Gawiser & Silk 
were able to find a model with G = 10 and low erg 

cluster abundance and 



which pro duced the correct z 
also that 



Verde et al. (2001) demonstrated that stronger 



constraints could be placed on G using the observed cluster 
X-ray size-temperature relation). 

Because of the prohibitively large amount of time re- 
quired for a full analysis of the likelihood contours for model 
parameters obtainable from SZ surveys (even with only four 
parameters allowed to vary) , we will estimat e likelihood con- 



tours using a Fisher-matrix analysis (e.g. Jungman et al 



Fo: the purposes of these examples, we consider only the 1996). This works well for the PLANCK survey, where the 



three cosmological parameters fib, fio and erg to be free (and 
assume fio -I- Ao = 1). The effects of varying h may be sim- 
ply scaled out of any results. These parameters are already 
well-measured from analysis of CMB experiments (e.g. |Net 



terfield et al. 2001 



) and cluster abundances (e.g. 



Eke, Cole 



& Frenk 1996). We will therefore adopt prior probabilities 



for the values of these parameters in our likelihood anal- 
ysis. Each parameter is assumed to have a Gaussian prior 
probability distribution with a width as given in Table ^ 
which also lists the true value of the parameter assumed 
in our likelihood analysis (we choose representative values. 



Gaussian assumption made in the Fisher-matrix analysis is 
justified. For the BOLOCAM survey however, this assumption 
is rather poor. To demonstrate this we show in Fig. ^ the 
likelihood contours obtainable for a 1 square degree BOLO- 
CAM survey of |S'i43ghz| > ImJy clusters, if redshifts are 
measured for all clusters. The thin lines indicate the results 
from the Fisher-matrix analysis (dashed lines assume no pri- 
ors, solid lines assume the priors from Table |^, marginal- 
ized over fio and erg. To compare to this, we generated a 
large, four-dimensional grid of models and then generated 
a large number of Monte-Carlo realizations of the observa- 
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Figure 7. Properties of SZ clusters for different sets of cosmological parameters. Dotted and daslied lines show results for rCDM and 
ACDM respectively. The left-hand panel shows counts as a function of flux at 143GHz (error bars indicate the expected cosmic variance 
in 1 square degree fields for the rCDM and ACDM cosmologies), the centre panel shows the redshift distribution of clusters brighter 
than |Si43GHzl = l™Jyi a-nd the right hand panel shows the distribution of angular sizes for the same clusters. All results are for the 
parameters of the BOLOCAM experiment and the J2000 P{y)- 




Figure 8. Properties of SZ clusters for different values of ag as indicated in the figure. The left-hand panel shows counts as a function 
of flux at 143GHz, the centre panel shows the redshift distribution of clusters brighter than |S'i43GHz| = ImJyi and the right hand panel 
shows the distribution of angular sizes for the same clusters. All results are for the parameters of the BOLOCAM experiment and the J2000 

Piv)- 



tional dataset. The best fitting model to the mock dataset 
was determined using a maximum likelihood approach, and 
the results used to construct likelihood contours in the G- 
Jlb plane. It is clear that the Gaussian assumption used in 
the Fisher-matrix analysis is a poor approximation, with 
the true likelihood contours depending on the parameters 
in a more complex way. Nevertheless, the Fisher-matrix ap- 
proach gives at least a rough idea of the constraints that 
can be placed on model parameters (and is very rapid to 
compute). 

Figure |ll| shows the constraints obtainable on the non- 
Gaussianity parameter G by the BOLOCAM and PLANCK ex- 
periments. In each panel we marginalize over two of the free 
parameters and show the error contours in the plane of the 
remaining two parameters. Thick and thin lines show re- 
sults with and without priors (priors are included only for 
the BOLOCAM experiment — for the PLANCK experiment the 
priors of Table |3| do not significantly alter the confidence re- 
gions). Dashed lines indicate results when no redshift infor- 



mation is available, while solid lines show results if redshifts 
are measured for all clusters. 

It is immediately apparent that, for the BOLOCAM ex- 
periment, no interesting constraints on G can be obtained 
unless either priors are assumed for the other parameters, or 
if redshifts are measured for all clusters (preferably both). 
(Note that in the left-hand panels the confidence ellipses for 
the case of no redshift information and no priors is so large 
that they he entirely outside the regions plotted.) With full 
redshift information and with the priors of Table ^ we find 
that BOLOCAM will measure G to an accuracy of ±40% (Icr). 
(It must be kept in mind, however, that the Fisher matrix 
analysis is only approximate for this survey — see Fig. p^ ) 
Interesting constraints on the three cosmological parame- 
ters considered can also be obtained if redshifts are available 
for all clusters in the bolocam sample, although there are 
strong degeneracies between the parameters (e.g. between 
(jg and Jlo in particular). The figure also demonstrates de- 
generacies between the various parameters, in particular be- 
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Figure 9. Properties of SZ clusters in ACDM for different degrees of non-Gaussianity. Solid, dotted and dashed lines show results for 
the RGS P{y) using G = 1 (i.e. Gaussian), 2 and 5 respectively. The left-hand panel shows counts as a function of flux at 143GHz, the 
centre panel shows the redshift distribution of clusters brighter than |S'i43GHzI = ImJy; and the right hand panel shows the distribution 
of angular sizes for the same clusters. All results are for the parameters of the BOLOOAM experiment and the ACDM cosmology. 



1^ 




-1 



Figure 10. Constraints in the G-f2b plane obtainable from the 
BOLOCAM survey assuming redshifts can be obtained for all clus- 
ters brighter than |S'i43GHz| = Imjy. Thin lines are the esti- 
mates from a Fisher matrix analysis, while thick lines were ob- 
tained from a large number of Monte-Carlo realizations of the 
dataset. Dashed lines show the results with no priors, while solid 
lines show results using the priors of Table ^ Contours show the 
68.3%, 95.4% and 99.7% confidence regions. The cross shows the 
true values of the parameters. 



tween ag and G. Increasing as or G tends to increase the 
total number of clusters of all fluxes and also to strongly 
increase the numbers seen at high redshifts (see Figs. ^ and 
^ leading to the degeneracy seen in Fig. 

For the PLANCK experiment simple number counts of SZ 
clusters will constrain G to an accuracy of ±5% while full 
redshift information would reduce that uncertainty to ap- 
proximately ±2% (since it is unlikely that redshifts would 
be available for all clusters the actual constraint on G will 



lie somewhere in between these two values). Even with- 
out redshift information, the PLANCK experiment will pro- 
vide extremely tight constraints on cosmological parameters 
through the SZ effect (which will supplement the determina- 
tions from the primordial CMB anisotropies which are the 
main goal of this experiment). It should be noted that, in 
the case of PLANCK the tight constraints shown here do not 
account for the sy stem atic uncertainties due to theoretical 
uncertainties (see §3.2), which would at present be the dom- 
inant source of error. 



3.3.3 Constraints on Preheating 

PreheatiiiE of clusters will affect the total number counts 
(see Fig. O where we show results for a range of value s for 
the minimum entropy parameter K34, defined in § 2.2.2 ), by 
dramatically reducing the number of faint clusters, while 
having very little effect on the redshift distribution of SZ 
sources. Preheating also increases the sizes of observed clus- 
ters, producing a tail to large 6 in the distribution of angular 
sizes. The effect is small and would require good statistics 
(and so a significantly larger survey than the 1 square de- 
gree proposed for bolocam), it is a clear signal of entropy 
injection since no other factor we examined has a similar 
effect on sizes. 

Figure |l^ shows confidence regions for the entropy pa- 
rameter K34 (the format follows Fig. ^). Again we see that 
useful constraints on K34 from the BOLOCAM experiment 
can be obtained only if redshift information is available or if 
priors are assumed for the various cosmological parameters. 
Assuming this to be the case, we find that BOLOCAM can con- 
strain K34 to an accuracy of ±60%. The vastly larger area 
surveyed by PLANCK will allow it to constrain K34, to an 
accuracy of ±6% assuming full redshift information. Note 
again that the effects on SZ statistics of altering K34 are 
degenerate with changes in all three of the cosmological pa- 
rameters considered in our analysis. 
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Figure 11. Left-hand panels: Joint constraints on f2o and Q^, (top panel), Co and G (middle panel) and Qq and erg (lower panel) obtainable 
from a 1 square degree BOLOCAM survey to |Si43GHzI = Imjy using a Fisher matrix analysis. In each panel we have marginalized over the 
two parameters not shown. Thin lines show results with no priors, while thick lines use the priors from Table ^ Dashed lines make use of 
only the number counts of SZ clusters, while solid lines make use of the their redshift distribution also (and assume that all clusters have 
measured redshifts). The contours show the 68.3% (i.e. Icr) confidence region. Right-hand panels: Same for a full-sky PLANCK survey to 
|5i43GHzl = Imjy. Note that we do not show lines including priors for the PLANCK experiment, since these do not significantly alter the 
confidence ellipses. 
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Figure 12. Properties of SZ clusters in ACDM for different levels of entropy injection. Solid, dotted and dashed lines show results for 
K34 = 0.0, 0.3 and 0.8 respectively. The dot-dashed lines showr the results when the entropy injected varies with the formation redshift 
according to K34 = 0.8/(1 + Zf). The left-hand panel shows counts as a function of flux at 143GHz, the centre panel shows the redshift 
distribution of clusters brighter than |5i43ghz| = Imjy, and the right hand panel shows the distribution of angular sizes for the same 
clusters. All results are for the parameters of the bolocam experiment and the ACDM cosmology with the J2000 P{y)- 



4 DISCUSSION 

We have made an in-depth analysis of the statistical prop- 
erties (including counts, redshift and size distributions) of 
clusters that should be detected through the thermal SZ ef- 
fect in forthcoming experiments, in particular BOLOCAM and 

PLANCK. 

Our model is based upon simple analytic calculations 
allowing us to explore a wide range of parameter space which 
is particularly important given that the details of the gas 
distribution within clusters arc unknown, and also allows 
for rapid calculation of confidence regions for cosmological 
and gas distribution parameters. Our standard calculation 
consists of a distribution of dark-matter halos with a mass 
function matched to that found in numerical simulations 
and with a distribution of formation redshifts motivated by 
the Press-Schechter theory. These halos are filled with gas 
in hydrostatic equilibrium and the Compton y parameter is 
calculated as a function of cluster radius. Finally, we account 
for smearing of the cluster SZ profile by the experimental 
beam and then compute the observable SZ flux above the 
self-consistently determined background. The result is the 
distribution of SZ clusters as a function of total SZ flux, 
redshift and angular size. 

We have quantified the effects of several theoretical un- 
certainties on the statistical quantities measurable from SZ 
surveys. We find that the uncertainty in how cluster forma- 
tion redshifts determine the thermodynamic properties of 
their gas leads to an uncertainty in cluster number counts 
comparable to the statistical uncertainty expected for the 
BOLOCAM survey. (Due to the broad redshift distributions 
characteristic of the SZ effect we find that sample variance 
is typically very small, i.e. close to the Poisson limit for the 1 
square degree BOLOCAM survey which we considered.) Other 
uncertainties, such as the exact form of P{y) and the density 
and temperature profiles of cluster gas have slightly smaller 
effects, but will still dominate over the statistical errors for a 
full-sky PLANCK survey. As such, it will be crucial to further 
our understanding of cluster gas properties before future SZ 
surveys can be exploited fully. 



Finally, we explored the constraints that may be set on 
cosmological parameters and also the presence of any non- 
Gaussianity in the initial conditions or preheating of ICM 
gas in the early Universe. Redshift distributions are poten- 
tially very sensitive to Qq, as and the presence of primor- 
dial non-Gaussianity. Preheating of clusters has much less 
effect on the redshift distributions but can significantly re- 
duce the number of faint SZ clusters and also produces a tail 
to large sizes in the SZ cluster size distribution. A Fisher- 
matrix analysis reveals that BOLOCAM will provide interest- 
ing constraints on cosmological parameters and the degree of 
non-Gaussianity (modulo uncertainties in K34) and/or pre- 
heating only if redshifts are measured for all the detected 
clusters (or if prior probability distributions are assumed 
for the cosmological parameters). The PLANCK experiment, 
due to the large area it will survey, should measure the 
non-Gaussianity or preheating parameters, G and K34 re- 
spectively, to accuracies of 5-10% even without redshifts for 
any of the clusters. Furthermore, the PLANCK SZ survey will 
provide independent and highly accurate constraints on cos- 
mological parameters that will complement those obtained 
from analysis of primordial CMB anisotropies. 

In conclusion, surveys in the SZ effect will soon provide 
valuable measurements of various cosmological parameters 
which will compliment measurements from other techniques. 
Assuming that the theoretical uncertainties highlighted in 
this paper can be adequately resolved, the SZ effect also has 
the potential to constrain the presence of non-Gaussianity 
and also to provide strong constraints on the distribution 
of gas within clusters and the presence or otherwise of any 
preheating. 
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Figure 13. Left-hand panels: Joint constraints on K34 and (top panel), K34 and Qq (middle panel) and K34 and erg (lower panel) 
obtainable from 1 square degree BOLOCAM survey to |S'i43ghzI = Imjy using a Fisher matrix analysis. In each panel we have marginalized 
over the two parameters not shown. Thin lines show results with no priors, while thick lines use the priors from Table ^. Dashed lines 
make use of only the number counts of SZ clusters, while solid lines make use of the their redshift distribution also (and assume that 
all clusters have measured redshifts). Right-hand panels: Same for a full-sky PLANCK survey to |Si43GHz| = Imjy. Note that we do not 
show lines including priors for the PLANCK experiment, since these do not significantly alter the confidence ellipses. 
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